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Abstract — In epidemic modeling, the term infection strength 
indicates the ratio of infection rate and cure rate. If the 
infection strength is higher than a certain threshold - which we 
define as the epidemic threshold - then the epidemic spreads 
through the population and persists in the long run. For a 
single generic graph representing the contact network of the 
population under consideration, the epidemic threshold turns 
out to be equal to the inverse of the spectral radius of the 
contact graph. However, in a real world scenario it is not 
possible to isolate a population completely: there is always some 
interconnection with another network, which partially overlaps 
with the contact network. Results for epidemic threshold in 
interconnected networks are limited to homogeneous mixing 
populations and degree distribution arguments. In this paper, we 
adopt a spectral approach. We show how the epidemic threshold 
in a given network changes as a result of being coupled with 
another network with fixed infection strength. In our model, 
the contact network and the interconnections are generic. Using 
bifurcation theory and algebraic graph theory, we rigorously 
derive the epidemic threshold in interconnected networks. These 
results have implications for the broad field of epidemic modeling 
and control. Our analytical results are supported by numerical 
simulations. 



I. Introduction 

In the existing individual-based epidemic models, the in- 
teraction and consequently the infection spreading process 
is driven by a single graph, the contact network by which 
individuals are in physical contact. However, in order to 
study epidemics in cyber-physical systems, a more elaborate 
description of the interaction is required. Several researchers 
from computer science, communication, networking, and con- 
trol communities are working on describing this complex 
interaction by using multiple interconnected networks |fl] , . 
The study of the spreading of epidemics in interconnected 
networks is a major challenge of complex networks, which 
has recently attracted substantial attention [3|-[6|. 

The N-Intertwined Mean-Field Approximated (NIMFA) 
model, first proposed by Van Mieghem Q, pointed out the 
specific role of a general network on the spreading process. 
Before NIMFA, most network-based epidemic models con- 
sidered aggregated networks characterized by a given node 
degree distribution Jg]. The NIMFA model has triggered 

F. D. Sahneh & C. Scoglio are with the Department of electrical and 
Computer Engineering at Kansas State University, Manhattan, KS, 66502. 
-E-mail: {faryad,caterina} @ ksu.edu 

F. N. Chowdhury is with the Directorate for Social, Behavioral & Eco- 
nomic Sciences, National Science Foundation, Arlington, VA, USA-E-mail: 
fchowdhu @ nsf . gov 



a pervasive amount of research on epidemic spreading on 
general networks, in different scenarios and with different 
compartments [9|-| 14 1. The key aspect of this class of models 
relay on the use of rigorous spectral analysis to determine the 
evolution of the epidemic. In particular, it is found that if the 
infection strength is higher than a certain threshold - which 
is defined as the epidemic threshold [7|, then the epidemic 
spreads through the population and persists in the long run. 
For a single generic graph representing the contact network 
of the population under consideration, the epidemic threshold 
turns out to be equal to the inverse of the spectral radius of the 
contact graph. The epidemic threshold provides a measure of 
the network robustness with respect to epidemics: the larger 
the epidemic threshold is, the more robust the network is, since 
more instances of infection strength will not spread in the long 
run. 

Current research efforts are directed to establish the impact 
of interconnected networks in spreading processes. Intercon- 
nection of networks can only make the system more vulnerable 
to infection propagation. Therefore, it is expected that the 
epidemic threshold of a network is not increased when it is 
connected to another network. Results for epidemic threshold 
in interconnected networks are limited to homogeneous mixing 
populations and degree distribution arguments. In two 
networks following the standard configuration model and 
interconnected with their own intranetwork are studied. Two 
possible scenarios are considered: strongly-coupled networks 
and weakly coupled networks. In strongly-coupled epidemics, 
either the epidemic invades both networks or not spread 
at all. In contrast, in weakly-coupled network systems, an 
intermediate scenario can happen where an epidemic spreads 
in one network but does not invade the coupled network. 

The objective of this paper is to study how much the 
interconnection can affect the robustness of the network when 
considering two general networks with network topologies 
expressed by their adjacency matrices. In this paper, we study 
the spreading process of a susceptible-infected-susceptible 
(SIS) type epidemic model in an interconnected network of 
two general graphs. First, we prove that interconnection always 
increases the probability of infection. Second, we find that the 
epidemic threshold for a network interconnected to another 
network with a given infection strengths rigorously derived as 
the spectral radius of a new matrix which accounts for the 
two networks and their interconnection links. We make use of 
algebraic graph theory to demonstrate our findings. 

The main contribution of this paper is the use of spectral 



2 



analysis to analyze epidemic spreading in interconnected net- 
works. To the best of our knowledge, this is the first time such 
approach has been used. As a result of our analysis on two 
general interconnected networks, we show that assumptions 
on the level of connectedness are not necessary to determine 
the epidemic threshold, and the evolution of the spreading 
process. Consequently, our results are rigorous and general, 
and reproduce results of [4 |, as a specific case. 

This rest of the paper is organized as follows. Preliminary 
tools in graph theory and a background in epidemic modeling 
is the subject of Section [TT] In Section Hill SIS epidemic 
spreading in two interconnected network is modeled. Main 
results on the epidemic threshold are provided in Section [IV] 
Finally, simulation results are available in Section [V] 

II. Preliminary and Background 

A. Individual-Based Epidemic Models 

Epidemic modeling has a rich history. Biological epidemi- 
ology has produced significant number of deterministic and 
stochastic models. These models have been successful in 
providing insights and deep understanding of the epidemic 
process phenomenon leading to successful conclusions about 
prevention and prediction of epidemics. In JT3], a stochastic 
epidemic model was studied for a well-mixed homogenous 
population. However, this assumption on the population ap- 
peared to be too simplistic in order to capture realistic cases. 
The theory of random networks was employed to generate 
models to represent contact patterns among individuals within 
a population. Specifically, results were reported in ifl&l for 
heterogeneous networks and in [8] for scale free networks. In 
the search for detailed and general models, individual-based 
epidemic models were proposed, where the contact network is 
represented by a graph. In particular, in the NIMFA model [7 1 , 
the probability of infection for each individual are the system 
states. For this model, the epidemic threshold is shown to be 
equal to the inverse of the spectral radius of the contact graph. 

B. Graph Theory 

Graph theory (see lUTl . 0~8]) is widely used for representing 
the contact topology in an epidemic network. Let Q — {V, £ } 
represent a directed graph, and V = {1, N} denote the set 
of vertices. Every agent is represented by a vertex. The set of 
edges is denoted by £ c V x V. An edge is an ordered pair 
G £ if agent i can potentially be directly infected by 
agent j. Mi = {j G V \(i,j) G £} denotes the neighborhood 
set of vertex i. Graph Q is said to be undirected if for any edge 
(i, j) G £, edge (j, i) G £ . In this paper, we assume that there 
is no self loop in the graph, i.e., (i, i) ^ £ , and the contact 
graph is undirected. A path is referred by the sequence of its 
vertices. A path V of length k between vq, Vk is the ordered 
sequence (vo,...,Vk) where (uj_i,i;j) 6 £ for i = 1, k. 
Graph Q is connected if any two vertices are connected with 
a path in Q. A = [ay] G M. NxN denotes the adjacency matrix 
of Q, where ay = 1 if and only if G £ else ay = 0. 
A graph is connected iff its associated adjacency matrix is 
irreducible. The largest eigenvalue of the adjacency matrix A 
is called spectral radius of A and is denoted by \i(A). 



A network of two interconnected graphs Q\ and Qi is 
represented by the set of non-overlapping vertices {V m }^ =1 
and the set of edges {£ mn }m,n=i> where £„ m CV m xV„ 
denotes the connection between vertices of V m to vertices of 
V n for each m,n G {1,2}. Connections between vertices of 
V m to vertices of V„ can be represented by matrices A mn , for 
each m, n G {1, 2}. 

In order to account for the hops of a path between the two 
interconnected graphs, we define the following class of paths. 
Without loss of generality, we assume that each path starts 
from i G Vi. 

Definition 1: A path from node % G Q% to node j is of class 
(Zi, l s ), with non-negative integers l\, l s , if it first make 
1% jumps in Q\ then goes to Q2 and make I2 jumps in Q2 then 
goes back to Q\ and makes I3 jumps in Q\ and so on until it 
makes the last l s jumps to reach j. 

It can be inferred from the above definition that a path of 
class (h, l s ), has length L = (s — 1)+ h + ■■■ + l s . 

III. Modeling SIS Spreading in Interconnected 
Networks 

In this paper, we study the spreading process of a 
susceptible-infected-susceptible (SIS) type epidemic model in 
an interconnected network of two graphs. In order to develop 
the model for the case of interconnected network, first we 
review the model for SIS spreading over a single graph. 

A. SIS Epidemic Spreading over a Single Graph 

Consider a network of N agents where the contact is 
determined by the adjacency matrix A. Agent j is a neighbor 
of i, denoted by j G Mi, if it can contract the infection to agent 
i. If j is a neighbor of i then ay = 1, otherwise ay = 0. In the 
SIS model, the state X.- L (t) of an agent i at time t is a Bernoulli 
random variable, where Xi(t) — if agent i is susceptible 
and Xiit) = 1 if it is infected. The curing process for 
infected agent i is a Poisson process with curing rate 5 G M + . 
The infection process for susceptible agent i in contact with 
infected agent j ' ^ i is a Poisson process with infection rate 
f3 G K + . The competing infection processes are independent. 
Therefore, a susceptible agent effectively becomes infected 
with rate j3Yi(t), where Yi(t) = J2jLi a ijXj(t) is the number 
of infected neighbors of agent i at time t. The ratio of the 
infection rate (3 over the curing i5 is the infection strength 
t = j. A schematic of SIS epidemic spreading model over a 
graph is shown in Fig. Q] 

Denote the infection probability of the i-th agent by pi = 
Pr[X;(t) = 1]. It has been shown that the marginal probabili- 
ties pi do not form a closed system. Actually, the exact Markov 
set of differential equations has 2 W states. Van Mieghem et. al. 
1171 used a first order mean-field type approximation to develop 
the NIMFA model, a set of ordinary differential equations 

N 

Pi = - Pi) ^2 a% '^i ~ Spi > i G N }' w 

3=1 

which represents the time evolution of the infection probability 
for each agent. There is no approximation on the network 



Fig. 1 . Schematics of a contact network along with the agent-level stochastic 
transition diagram for agent i according to the SIS epidemic spreading model. 
The parameters and <5 denote the infection rate and curing rate, respectively. 
Yj (t) is the number of the neighbors of agent i that are infected at time t. 




topology in this model. According to this model, it is proved 
that if the infection strength r = (3/5 is less than the 

initial infection probabilities die 



threshold value t c = 
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out exponentially. If the infection strength r = (3/8 is higher 
than t c , then infection probabilities will go to non-zero steady 
state values. 

B. SIS Epidemic Spreading over Interconnected Networks 

Consider two groups of agents of sizes Ni and N%. In order 
to facilitate the subsequent developments, we label the agents 
of the first graph Q\ from 1 to iVi, and the agents of the second 
graph Qi from Ni + 1 to N\ + N2. The collective adjacency 
matrix A, defined as 



A 



An 
A 2 i 



Al2 
A 2 2 



(2) 



represents the contact between all of the agents. Since the 
contact topology in this paper is undirected, An and A22 are 
symmetric matrices and A21 — A\ 2 - According to definition 
©, agent i is connected to agent j iff (A)ij = 1. A 
schematic of the interconnected contact network of the agents 
is represented in Fig. [2] 

The SIS spreading model over a single graph described in 
Section IIII-AI can be generalized in the following way. The 
curing rate for agents of graphs Q\ and Q2 are <5i S K + and 
82 <E K + , respectively. The infection rates /?n, /?i2, /?2i, ^22 £ 
R + are such that a susceptible agent of graph Q m receives 
the infection from an infected agent in Q n with the infection 
rate (3 mn , for m,n £ {1,2}. Similar to (Q~|l, the infection 
probabilities of the agents evolve according to the following 
set of differential equations: 



N t +N 2 



Pi = (1 - Pi){Pu ^2 aijpj + P12 ^2 a,,/',} - 5iPi, 

ie{l,-^i}, (3) 

Pi = (1 - Pi) {An ^2 ai 3Pl + X! ai 3Po} - S 2Pj, 

»6{JV 1 + l,...,JVi+J\T 2 }. (4) 



Fig. 2. A schematic of the coupling between graphs Q\ (in black) and Q2 
(in red), which each represent a contact network. The blue links represent the 
coupling between the nodes of the two graphs. Q\ and Q2 are not necessarily 
connected. However, the whole interconnected network is connected. 



Remark 1: Infection process is the result of interaction 
between a pair of agents. Therefore, actually the infection 
rate in (Q} equals to {3 = [lit where /1 6 R + is the rate that 
an infected agents transmits the infection and n 6 [0, 1] is 
the probability that a susceptible agent receives a transmitted 
infection. Similar arguments show that the four infection 
rates fin, (3i 2 , /J21, (322 in ® and (0]l are not completely 
independent of each other. Having (3n = Mi^i and ^22 = 
H2K2, the infection rates (3\2 and /?2i will have a form of 
(3\2 = a^i7T2,/?2i = Oifj,2^i, where a g R + is a positive 
scalar accounting for heterogeneity of interconnection and 
interaconnection. Therefore, the following constraint exists 
among the infection rates 



/3n/?22 = « 2 /3i2/?: 



21- 



(5) 



Comparing (Q~|i and (0, it can be concluded that intercon- 
nection increases the probability of infection. This conclusion 
is actually intuitive: when interconnected with other agents, 
there is more possibility to receive the infection. 

IV. Main Results 

A. Problem Statement 

Suppose that agents of graph Q\ are connected to agents of 
graph Q2, and the overall contact among the agents is deter- 
mined by A defined in (fJJ, where the following assumption 
holds for T22. 

Assumption 1: If there is no interconnection, infection can- 
not survive in Q2, i.e., 

1 

(6) 



>22 



S 2 K Ai(A 2 2) 

Under Assumption [TJ for the infection strength t\\ =0, 
the steady state value of infection probabilities of <[3j and 
© are necessarily zero for each agent. In this paper, we 
find a threshold value t u c such that for infection strength 
T11 > Tn, c the steady state infection probabilities take positive 
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values. Since interconnection always increases the chance of 
receiving the infection, we expect ru c < l/Ai(An). 

B. Equation for Epidemic Threshold 

We use bifurcation theory to find the epidemic threshold. 
From (0 an d the equilibrium points of the infection 
probabilities satisfy the following set of algebraic equations 



Pi 



Ni 



Nt+N 2 



Tn ^2 ai lP*3 +Tl2 ^2 ai lP*3 



J=A'i+l 



»G{l,...,JVi}, (7) 



N\ Nx+N 2 
T21 ^ ai l'P*3 + T22 ^'P/ 



J=JVi+l 

%e{Ni + l,...,Ni+N 2 }, (8) 



where, 



A A Pl2 A /?21 A 022 

Til — -T—,Tl2 — -T— ,T21 — —£—-,T22 — -J— ■ 
Ol Ol 02 C>2 



(9) 



Lemma 1: If the overall contact network is connected, the 
steady state values of the infection probabilities are either zero 
for all of the agents or absolutely positive for each agent. 

Proof: The idea of the proof is inspired from [7|. The 
steady state values for the infection satisfies and (0. 
Therefore, p* = for Vi G {1, 2Vi + iV 2 } is a solution for 
the steady state infection probabilities. Suppose there exists a 
node j such that p* > 0. According to (0 and (0, for any 
node i that is a neighbor of node j, i.e., dij / 0, the steady 
state infection probability is 



1 + Tn 2^=1 a ijPj + T 12 2^ 3 =N 1 + 1 a i]Pj 

if i G {1, ...,iVi} and 

* _ T 21 Ej=l aijPj + T 22 Lj=N 1+ l a ijPj 



(10) 



" 1 + T2X ESx «iiPj + ^22 E^r+! ^ ' 

if i G {ZVi + l,...,iVi + A^}, which is positive because 

2~2jli a i]Pj > or Ysf^N^+i a ijPj > °- Same procedure 
can be applied to the neighbors of node i, and so on. Hence, 
if the contact network is connected and at least one of the 
agents have nonzero infection probability, then p* > for all 

ie{l,...,Ni + N 2 }. ' ■ 

Before the epidemic threshold, origin is the only solution 
to (0 and (01. Epidemic threshold is the critical value m )C 
such that a second equilibrium point starts leaving the origin. 
A corollary of Lemma0is that the epidemic threshold t\\ c is 
such that p* = and J^L > for every i G {1, Ni + N 2 }. 
Taking the derivative of (0 and (0 with respect to Tn at 
rn = Tn )C and p* = yields 



N x +N 2 



M 

dm 



Ni 



'9tij 



N!+N 2 
+ T12 ^ a 



13 



9m ' 



3=JVi+l 

ie{i,...,JVi}, (li) 



9ru 



3=1 



T22 



9rii' 



ie{Ni + l,...,Ni + N 2 }. (12) 



Defining Vi = [J^ 



dp 



»i IT 



and V 2 = 

[^±i,..., 9p ^+ JV2 ] r , the equations QT) and fl2) can 
be equivalently expressed in the collective form as 



(13) 



The critical value of the infection strengths are those for which 
the above equation has a positive solution. Equation dT3l > can 
be written as 



Tll.c^ll T12A12 




'V{ 




vr 


T 2 lAj 2 T 2 2^22_ 




Y2. 




v 2 



Tii, c AnVi + T12A12V2 = Vi, 

T 2 lAj 2 Vi + T22A22V2 = V 2 . 



(14) 
(15) 



According to Assumption0 if Vi is positive then V2 = T21 (I— 
T22^-22) _1 Aj 2 Vi exists and is non-negative. Therefore, (fT3l is 
equivalently expressed as 



i?Vi = Vi 



where H is defined as 



# - Tii.cAn + t 21 t 12 A 12 (I - T22A22) Y A\ 



(16) 



(17) 



C. Effect of Coupling on Epidemic Threshold 

The rest of the analysis is to find the threshold value m c 
such that ([Tol l has a positive solution for Vi. The following 
results facilitate the proof of Theorem which is the main 
result in this paper. 

Lemma 2: The number of paths of length L from node i G 
Qi to node j corresponding to the class (Zi, ...,l s ) is: 

« the (i,j)-th entry of A l ^ l Ai2A%A 2 i ■■■A 2 - L A l {, if j G 
{l,...,iVi}, 

* the (i, j - 7Yi)-th entry of A l 1 1 1 Ai2A^ 2 A 2 i ■ • • Ai 2 A l 2 s , if 
j G {JVi + l,...,iVi + Ar 2 }, 

where A\ x = In iX Ni an d ^22 = In 2 xN 2 > by convention. 
Proof: We use induction for the proof. For L = 1, the 
number of paths from node i to j is equal to 1 if i is connected 
to j, and is zero otherwise. If j G {1, ...,N±}, path of length 
L = 1 corresponds to the class (1). Therefore, the number of 
paths from node i to j is equal to the (i,j)-th entry of An. 
If j G {A^i + 1, Ni + N 2 }, then a path of length L = 1 
corresponds to the class either (0, 0). In this case, the number 
of paths from node i to j is equal to the — iVi)-th entry of 
A12 — AiAi 2 A 2 . Therefore for L = 1, the Lemma is correct. 

Assume that for L = Lq the lemma statement is correct. 
Consider the first case where j G {1, ...,N{\. A path of length 
L = Lq + 1 from i to j is either of the class (Ji, l s + 1) 
or (7i, l s , 0). Such a path can be constructed from paths of 
length L from i to k of the class (Zi, ...,Z S ) then connected 
to node j from node k. 
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If the path from i to j is of class {ly, l s + 1), then the 
number of such paths is 



Y,(A[\A 12 A l &A 21 ■ ■ ■ A 21 A l {) lk {A 1 ) kj = 



k=l 



A[\A 12 A[%A 21 ---A 21 A l { +1 . 



If the path from i to j is of class (h, l s ,0), then the number 
of such paths is 



N!+N 2 

E 

fc=jVi+l 



(A[\A 12 A[lA 21 ■ ■ ■ A 21 A l {) l(k _ Nl) {A 



12) (k-Nt)j 



= A[\A 12 A[%A 21 ■ ■ ■ A 21 A l { A 12 

= A[\A 12 A[ 2 2 A 21 ■ ■ ■ A 21 A[>A 12 A° 2 . 

Hence, the theorem statement is correct for L = L +l and j £ 
{1, Ni}. Similar procedure can be followed to conclude the 
same result for j 6 {iVi + 1, Ni + N 2 }. ■ 
Theorem 1: The matrix Ht defined as 

H T 4 A n + a 2 T 22 A 12 (I - T 22 A 22 y x A T 12 . (18) 

is irreducible if the coupled network is connected. 
Proof: We show that 

Na-l 

S T = A u + A 12 A\ 2 + A i2A k 22 A T l2 (19) 
fe=i 

is irreducible. If Ht is shown to be irreducible, then 
A n + ol 2 t 22 A 12 A\ 2 + a 2 r 22 T k 2 A 12 A k 2 A\ 2 is 

irreducible. And hence, Ht = An + o?t 22 A\ 2 A\ 2 + 
u 2 T 22 Y^ =1 T 2 k 2 A X2 A k 2 A\ 2 = An + a 2 r 22 A 12 (I 
t 22 A 22 )~ 1 A{ F 2 is irreducible and the proof is completed. 

If Q\ is a connected graph, then An and as consequence 
Ht is irreducible. Assume that An does not represent a 
connected graph. Therefore, there exists a pair i, j such 
that there is no path between them in Qy. However, since 
the whole interconnected network is connected, there ex- 
ists a path from i to j. Suppose, the path is of class 
(h,lj2,i,li,2,l2,2,:.;h, s ,h,s+i), i-e., it makes h t i jumps 
in Qi to reach vertex k° ut , then it leaves Qi and enters 
Q2 and makes l 2< \ jumps in Q2, then enters Q\ at vertex 
k™ . This process goes on until it makes li :S +i jumps in 
Qi from k" L to reach vertex j. Matrix Ht is proved to be 
irreducible if we show that (fc°"*,fc™)-th entry of Ht is 



„out 



to 



positive for u = 1, s. Since, there is path from k 
fc™ which is of the class (0,/ 2 ,„,0), the (k° u ut , fc™)-th entry 
of A\ 2 A l 22 A\ 2 > l,because it is the number of such paths 
according to Lemma |2] As a consequence, (k° ut , fc™)-th entry 
of Ht is positive and therefore Ht is irreducible. Hence, the 
proof is completed. ■ 
Theorem 2: For 

1 



n - " MiJrt • <20) 

where Ht is defined in ( fT8l . the equation (foi l has positive 
solution for V\ and V 2 . That is Tn.c is the epidemic threshold. 
Proof: According to (0 and the definitions ©, we have 



Substituting for t 2 \T\ 2 in ( TT7| >. equation ( fT6*b gets the form 

tii,cH t Vi = V u (22) 

where Ht is defined in ( fT8l . In order for d22t to have 
solutions, Tn, c must be the inverse of one of the eigenvalues of 
Ht- However, the corresponding eigenvector V\ must have all 
positive entries. According to TheoremQ] Ht is an irreducible 
matrix. Therefore, V\ is positive only for the eigenvector 
corresponding to the largest eigenvalue. Therefore, Tn.c is 
equal to l/\i(H T )- ■ 

D. Interconnection Topology and Epidemic Spreading Modes 

We used a bifurcation method to find an expression for 
the epidemic threshold. According to our definition, the epi- 
demic threshold is a critical value such that for any infection 
strength t\\ > Tn.c, the steady state values of the infection 
probabilities are positive. Consider the special case where 
the infection strength t 22 is very close to the spectral radius 
M(A 22 ), i.e., T22 -> Ai(A 2 2). According to (Q~8]l and d20l >. 
the epidemic threshold tu )C —> as far as the whole contact 
network is connected. This argument is true even for very weak 
interconnection between the two networks Q\ and C?2- The 
reason for this observation is that since t 22 — > Xi(A 22 ), only 
a small amount of interconnection will lead to an outbreak in 
Q 2 . Since, the probability of infection in Q 2 becomes a positive 
value, according to Lemma Q] the probability of infection in 
Qi is also positive. Therefore in the case of t 22 — > Xi(A 22 ) 
and weak interconnection, the positive infection probabilities 
in Qi for small values of t\\ is only due to epidemic outbreak 
in Q 2 . 

The numerical simulations in Section [V] illustrates three 
possible curves of r n c as a function of t 22 , as shown in 
Fig. [3] Here, the blue curve belongs to the case of weak 
interconnection between the two graphs. As can be seen, the 
decrease in the epidemic threshold Tn.c is very slow for small 
values of t 22 , while there is a quite sharp drop in the values of 
Tn.c as T22 — ) Xi(A 22 ). For strong interconnection topology, 
shown by the green curve, the value of Tn.c is decreasing 
quickly for small values of t 22 . In this case, the infection 
in Qi starts to grow not only as the result of receiving the 
infection from Q 2 , but also as the result of a survivable internal 
infection force. The red curve is an intermediate between the 
two spreading modes. 



Theorem 3: The derivative 



drxi 



dr-2-2 



at r 22 = is 



dm 



dT 2 



\AT 2 xi 



1"22=0 



(23) 



where xi is the eigenvector of An belonging to \i{An) 
Proof: The matrix Ht from ( fT8l can be written as 



Hi 



An 



a 2 T 22 A 12 A{ 2 



o(r 22 ). 



mm = a Tn. c T 22 . 



(21) 



Therefore, taking the derivative of (1221 with respect to t 22 at 
T22 = yields 

^iiH + . \ (a 2 A 12 Aj 2 ) Xl 
dr 2 2 M(An) 

M{An) uT 2 2 



Multiplying ( 1241 by xf from left, we get 



l,c T 



dr 22 



x x A xx x\ 



Ai(A n ; 



12 Al 2 )xi 



1 



-An -7) 



Ai(An) dr 22 

In (f25j), since x\is the normalized eigenvector of 
responding to A 1 (A 11 ), we have x x AxiXi = Ai(^ 

(177^773^11 ~ -0 = f° r ^11 i s symmetric. T 
equation (125b becomes 



Ai(A 



Hence, 



11 



dT 22 Ai(-Au) 



(^i) T (^i) 



is found to be d2~3T >. 
According to (l23l and the proceeding arguments, \ 
interconnection topology measure 



\A\ 2 x x 



Ai(An)Ai(A 22 ) 

If Q,{Qi,Q 2 ) < 1, then for the infection strength 
above the threshold r n c in ( f20l >. the positive infection prob- 
ability in Qi is mostly due to external infections from Q 2 . 
While if £l(Gi,G 2 ) > 1, then for the infection strength t\\ 
right above the threshold t\\ c in d20l i, the positive infection 
probability in Q\ is mostly due to a survivable internal infec- 
tion force as the result of increased effective level of contact 
among agents of Q\. 

Remark!: A very interesting property of 0(^1,^2) de- 
fined in d27b is that it is a purely topological measure and 
does not depend on the epidemic parameters. 

V. Numerical Simulation Results 

We have generated two graphs according to the small world 
random network model [19|. The first network has N\ — 500 
vertices with Watts and Strogatz parameters for mean degree 
K\ = 10 and the rewiring probability j3i = 0.2. For this 
graph, the spectral radius is found to be Xi(An) = 22.0586. 
The second network has A^ 2 = 100 vertices with the Watts and 
Strogatz parameters for mean degree K = 2 and the rewiring 
probability fi 2 = 0.1. For this graph, the spectral radius is 
found to be Ai(A 22 ) = 4.3. For the interconnection of these 
two graphs, we use the following rule. All the potential edges 
between the first graph and the second graph are active with 
some probability ui, to be chosen. 

In the first simulation, f c i = Ai(An)rn )C is plotted as a 
function of f 2 = Ai(A 22 )t 22 , for three different values of 
10 = 0.01, 0.042, 0.2. The numerical results for the three cases 
are shown in Fig. [3] As can be seen from the blue curve, 
which is for u) — 0.01, for weak interconnection among the 
graphs the emergence of positive steady state values for the 
infection probability in graph Q\ is due to an outbreak in Q 2 . 
While for strong interconnection u> = 0.2, which is shown with 
the green curve in Fig. [3] the emergence of positive steady 
state values for the infection probability in graph Q\ is more 
due to increased level of effective contact among the agents 
in Q\. The red curve in Fig. [3] belongs to an intermediate 
interconnection strength, here u> — 0.042. 
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Fig. 3. Normalized epidemic threshold r c i = Ai(j4h)th. c of graph Q\ as 
a function of the normalized infection strength -T2 = Ai(A22)t"22 of graph 
(?2- The interconnection in (a) uj = 0.01, the blue curve, (b) u) = 0.042, 
red curve, and (c) u) = 0.2, green curve. A %50 reduction of the epidmeic 
threshold is observed for the normalized infection strengths (a) T2 = 0.925, 
(b) f 2 = 0.5, (c) f 2 = 0.05. 



According to Fig. |3]a %50 reduction of the epidemic 
threshold is observed in Q\for (a) uj = 0.01 and f 2 = 0.925, 
(b) uj = 0.042 and f 2 = 0.5, (c) u = 0.2 and f 2 = 0.05. 
We have plotted the curves of f>* = jj~ Yli=i Pi as a 
function of r 11 Ai(A 11 ). We have found the equilibrium values 
of p* by solving the algebraic equations ((T) and ((H). The 
numerical method for solving these equations is presented in 
the Appendix. 

VI. Conclusion 

In epidemic modeling, the term infection strength indicates 
the ratio of infection rate and cure rate. If the infection strength 
is higher than a certain threshold then the epidemic spreads 
through the population and persists in the long run. For a 
single generic graph representing the contact network of the 
population under consideration, the epidemic threshold turns 
out to be equal to the inverse of the spectral radius of the 
contact graph. However, in a real world scenario it is not pos- 
sible to isolate a population completely: there is always some 
interconnection with another network, which partially overlaps 
with the contact network. We study the spreading process of 
a susceptible-infected-susceptible (SIS) type epidemic model 
in an interconnected network of two general graphs. First, we 
prove that interconnection always increases the probability of 
infection. Second, we find that the epidemic threshold for 
a network interconnected to another network with a given 
infection strengths rigorously derived as the spectral radius of 
a new matrix which accounts for the two networks and their 
interconnection links. The main contribution of this paper is 
the use of spectral analysis to analyze epidemic spreading in 



7 



0.7 




Ai(An)rii 



Fig. 4. The mean steady state infection probability in G± as a function of 
the normalized infection strength of Ai(An)rn for graph G\. Black curve 
corresponds to the case where there is no interconnection. In this case, the 
epidemic threshold is tji c = l/Ai(An). All the other curves correspond 
to the case where T\\, c = ^ X l/Ai(An). For (a) the blue curve ui = 0.01 
and T2 = 0.925, (b) the red curve ui = 0.042 and T2 = 0.5, (c) the green 
curve ui = 0.2 and T2 = 0.05. 

interconnected networks. Our results have implications for the 
broad field of epidemic modeling and control. 
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VIII. Appendix 

In order to numerically solve d7} and (|8), we can use the 
recursive formula 

Ni ,,s /Vi+JV 2 

ieft-,*}, (28) 

Ni ,,s /V1+JV2 

iG{iVi + l,...,JVi + JV 2 }, (29) 

which converges to the fixed points y* :— ^ , , for i G 
{l,...,Ni + A 2 }. The steady state values of the infection 
probabilities are then p* = ^ , . 



